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Introduction 

Because  of  continuing  current  interest,  the  properties  of  gases  and  gas  mixtures 
of  interest  in  chemical  laser  applications  were  studied.  This  final  technical  report 
constitutes  the  second  in  this  series  and  supplements  those  substances  and  mixtures 
reported  earlier  (1]  *.  The  present  report  also  recommends  further  work. 


Theoretical  and  Empirical  Methods 


In  the  earlier  report  the  approach  was  taken  that  the  Lennard-Jones  6-12 
potential  was  used  to  calculate  the  viscosity  and  diffusion  coefficients.  For  user  con¬ 
venience  a  brief  summary  of  that  material  is  reproduced  here.  For  gaseous  viscosity 
the  mixture  viscosity,  is  given  by 

1+  Z  ... 

^M"X+Y  (1) 


where 

xx2  2xt  Xj  Xj8 

X  = —  + - +  — 

Mil  Ml2  M22 

3Ai2  jxi2  Mj  2x1x2  (M,+ M2)2  Xj2 
Y  =  ~5~fe  4M,M2 

3A«  r  Mi  r(Mj  +  M2)2  m12  Mj-. 

Z=~  4M.M,  -1J+’^  m;/ 


In  these  equations  Mlf  M2,  xlt  and  are  the  molecular  weights  and  mole  fractions 
of  components  1  and  2  and  the  jj12,  and  ji22  represent  the  viscosities  of  component 
1,  a  hypothetical  gas  12  and  component  2.  The  gas  12  is  assumed  to  obey  the  same 
potential  as  gases  1  and  2  but  to  have  a  molecular  weight  M12  of  2Mt  M2/(  Mt  +  M2) , 
a  collision  diameter  aJ2  of  (Oi  +  a2)  /2  and  a  well  depth  €  12A  of  iA)  (€2/k).  A12 
is  a  collision  integral  function,  tabulated  simply  as  A  in  |1) .  The  are  given  by 

10'  =  266.  93^7/{o'  On  (kT/Cy)},  {} ;  };|}  (2) 

where  the  022  is  a  collision  integral,  described  and  also  partially  tabulated  in  ( 1] . 

The  quantities  tabulated  in  our  tables  are  /jtl,  fjtz,  and  ^22  so  that  equation  (1)  should 
then  be  used  to  calculate 


*  Numbers  in  square  brackets  refer  to  the  reference  list. 
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For  diffusion  the  equation 


101  D«  - 


262.84  gT  1  cm1 

Po*.^  V^ij  nij  (kT^V  sec 


(3) 


can  be  used  and  the  quantities  tabulated  are  Dn,  D12,  and  Dj,,  for  P  =  1  atm. 

For  thermal  conductivity  a  more  empirical  approach  was  adopted.  Values  for 
the  pure  components  were  taken  from  earlier  works  [2,3]  where  possible.  These  val¬ 
ues  basically  are  curve-fits  for  the  lower  temperatures  (i.e. ,  below  about  1000  K) 
faired  into  power  function  fits  to  the  data  (if  any  existed)  above  1000  K  or  from  a  sim¬ 
ilar  analysis  of  high  temperature  viscosity  data.  Values  for  the  mixtures  are  merely 
the  average  (i.e. ,  equimolar)  values  for  the  pure  components,  based  on  the  premise 
that  the  thermal  conductivity  of  a  mixture  is  a  linear  function  of  its  composition. 

The  summary  of  theoretical  work  in  our  earlier  report  led  to  the  conclusion 
that  no  intermolecular  potential  is  capable  of  representing  the  data  for  all  properties 
to  within  the  experimental  error.  This  is  due  to  orientation  effects  and  inelasticity 
in  collisions  coupled  with  the  inability  to  perform  the  mathematical  calculations  to 
sufficient  accuracy  and  the  fact  that  high  energy  collision  and  spectroscopic  data  may 
yield  still  different  intermolecular  potential  values.  This  conclusion  is  in  part  sup¬ 
ported  by  the  more  recent  analysis  [6,7]  of  the  ability  of  the  m-6-8  potential  to 
simultaneously  fit  viscosity  and  second  virial  coefficient  data  for  non-polar  polyatomic 
gases.  It  was  found  that  additional  terms  were  necessary  to  represent  the  virial  co¬ 
efficient  behaviour.  This  implied  that  a  more  comprehensive  potential  should  then 
likewise  be  used  to  calculate  the  viscosity  (and,  of  course,  other  properties).  This 
is  an  enormous  complication,  both  from  the  fundamental  difficulty  of  calculating  func¬ 
tions  for  comprehensive  potentials  and  from  the  more  practical  viewpoint  of  fitting 
such  potentials  to  the  experimental  data.  While  it  is  stated  [6]  that  "work  is  in  prog¬ 
ress  to  modify  the  m-6-8  potential  to  include  non-spherical  effects”  one  reads  on  the 
same  page  that  "one  assumes  that  the  molecules  collide  with  a  fixed  relative  orienta¬ 
tion  and  collision  integrals  are  calculated  for  each  orientation  independently.  The 
final  values  of  the  collision  integrals  are  then  determined  by  a  statistical  averaging 
over  all  these  possible  orientations".  Here  is  the  important  fact  that  the  calculations 
to  date  assume  a  fixed  relative  orientation.  This,  coupled  with  the  averaging  require¬ 
ment  will,  in  the  opinion  of  the  writer,  introduce  sufficient  uncertainty  so  as  to  negate 
possible  improvements  by  modifying  the  potential.  Indeed,  it  appears  probable  that 
this  stage  has  already  been  reached  from  a  fundamental  viewpoint  by  the  failure  of 


the  potential  to  correctly  model  the  dispersion  coefficients,  as  was  already  mentioned 
in  the  discussion  of  equation  (4)  of  our  earlier  report  [1] .  An  extensive  review  (81 
of  some  aspects  of  intermolecular  forces  does  nothing  to  dispel  the  conclusion  that 
realistic  potentials  are  not  yet  available  for  molecules  of  complex  structure.  Hence, 
in  the  present  work  our  previous  use  of  the  Lennard-Jones  6-12  potential  is  continued, 
especially  as  no  extensive  investigations  of  the  ability  of  other  potentials  to  represent 
the  transport  properties  of  gaseous  mixtures  appear  to  have  been  made. 

Selection  of  Lennard-Jones  6-12  Parameters 

The  parameters  for  substances  which  have  previously  received  attention  in 
the  literature  were  taken,  where  possible,  from  the  tables  of  Svehla  (10)  for  reasons 
detailed  in  our  earlier  report.  Parameters  for  other  substances  were  estimated  using 
procedures  either  identical  or  similar  to  those  of  Svehla.  The  following  table  repro¬ 
duces  values  of  parameters  for  pure  component  substances  not  tabulated  in  our  previous 
report. 

Lennard-Jones  6-12  Parameters  for  Pure  Gases 


Gas 

e/k<°K) 

a  (1) 

i/Vm 

He 

10.22 

2.551 

0.49980 

N, 

113.50 

3. 566 

0. 18894 

A 

138.20 

3.287 

0. 15822 

59.70 

2.827 

0. 70430 

D 

35.20 

2.664 

0. 70712 

NH, 

558. 30 

2.900 

0. 24232 

NDj 

337.00 

2.900 

0. 22357 

n2h« 

760.00 

3.180 

0. 17665 

n2d4 

462.00 

3.380 

0.16664 

*2 

112.60 

3.357 

0. 16223 

F 

112.60 

2.968 

0. 22942 

NFS 

175. 00 

4.154 

0. 11868 

n2f4 

240.00 

3.880 

0. 09805 

CIFj 

448. 30 

4.900 

0. 08770 

HF 

330.00 

3.148 

0. 22357 

DF 

199.10 

2.826 

0.21223 

HC1 

344.70 

3.339 

0.16559 

DC1 

208.00 

2.980 

0. 16340 

H 

37.00 

2.701 

0.99210 

Dj 

35.21 

2.952 

0. 50000 

No  analogous  table  to  that  of  our  previous  report  for  all  the  gas  pairs  considered  here 
is  given  in  this  report  as  the  tables  of  the  present  report  were  generated  from  a  com¬ 
puter  program  in  which  the  mixture  parameters  were  internally  generated  from  the 
pure  co’  oonent  vab  ts. 
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Calculation  of  Transport  Properties  of  Purv  Gases  and  of  Gas  Mixtures 

The  computer  print-out  was,  in  fact,  programmed  so  that  the  tables  of  viscosity 
and  diffusion  coefficient  appearing  in  this  report  are  edited  copies  of  the  computer 
sheets.  The  thermal  conductivity  values  were  based  on  experimental  data,  correlated 
with  theory  if  possible  as  described  in  [3] ,  or  were  generated  by  applying  corrections 
for  internal  degrees  of  freedom  to  the  calculated  translational  effect.  In  the  latter 
case,  some  substances  showed  a  significant  divergence  between  the  theoretical  and 
the  experimental  values.  In  such  cases,  the  calculated  values  were  reduced  to  bring 
them  into  general  accord  with  the  experimental  data.  A  comment  is  made  on  this  in 
the  Recommendations  section.  In  our  programming  we  adopted  the  fittings  of  Neufeld, 
Janzen,  and  Aziz  [11]  to  generate  the  6-12  collision  integrals.  It  should  be  noted  that 
their  representation  is  only  valid  for  reduced  temperatures  from  0.  3  to  100  and  not 
from  0. 1  to  0. 3  or  from  100  to  400.  Some  problems  were  encountered  before  this 
fact  was  realized.  While  the  lower  end  gives  no  practical  difficulty,  substances  with 
low  characteristic  well  depths,  €/k ,  can  yield  reduced  temperatures  over  100.  A  suit¬ 
able  blocking  procedure  was  thus  written  into  the  program  to  prevent  the  computer 
extrapolating  the  calculation  into  regions  where  the  reduced  temperature  exceeded 
one  hundred.  For  the  few  cases  where  this  occurred,  tables  were  numerically  gen¬ 
erated  using  the  approximations 


n11=  1.10670  (t*)-°-15?s 
022  =  0.07183  (T*)-°-t62S 


(4) 


which  have  been  found  to  yield  collision  integrals  accurate  to  1/4  percent  for  reduced 
temperatures  over  25.  In  fact,  the  use  of  the  Lennard-Jones  potential  at  such  high 
reduced  temperatures  is  physically  unrealistic  due  to  the  predominance  of  the  repul¬ 
sive  part  of  the  potential  for  which  an  exponential  dependence  on  intermolecular 
separation  is  preferable  to  a  power  one. 


Summary  and  Recommendations 

In  the  eighteen  months  duration  of  this  project  the  first  six  months  saw  the 
preparation  of  tables  of  the  thermal  conductivity,  viscosity,  and  Ficktan  diffusion  co¬ 
efficient  for  38  binary  mixtures  and  9  pure  substances  which  were  distributed  through 
AFOSR  [1] .  The  present  report  gives  tables  for  an  additional  69  binary  mixtures. 
Both  sets  of  tables  refer  to  the  dilute  gas  state  and  present  tables  at  100  K  increments 
from  100  K  to  3000  K.  The  division  of  work  between  the  two  sets  of  tables  followed 
AFOSR  assessment  of  relative  priority  of  the  different  mixtures  and  pure  gases. 
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As  a  result  of  the  work,  some  facts  and  conclusions  can  be  stated.  These  can 
be  summarized  as 

1.  No  tables  have  yet  been  generated  for  pressures  above  the  atmospheric 
pressure  (=  dilute  gas)  state.  While  a  detailed  analysis  would  require  a 
prolonged  effort,  an  approximate  study  could  be  made  using  empirical  tools 
of  moderate  accuracy.  Some  of  the  problems  here  are  connected  with  the 
accuracy  with  which  AFOSR/AFWL/others  might  require  the  effect  of  pres¬ 
sure  to  be  known. 

2.  Of  the  three  properties  with  which  we  have  been  concerned,  thermal  con¬ 
ductivity  is  the  most  uncertain.  While  estimates  of  the  intermolecular 
potential  affect  all  three  properties,  for  thermal  conductivity  above  the 
effect  of  the  internal  degrees  of  freedom  can  outweigh  the  translational 
energy  contribution.  It  is  proposed  to  pursue  the  problem  of  estimating 
the  internal  effect  by  generalized  correlations  firstly  and  then  by  rigorous 
theory  if  time  allows. 

3.  Only  one  intermolecular  potential  function  has  been  used  to  any  extent  to 
calculate  mixture  propertied.  Attention  should  be  expanded  to  other  potentials. 

4.  A  useful  further  approach  would  be  to  use  the  Prandtl  number  as  a  means 
of  obtaining  thermal  conductivies  based  on  critically  evaluated  viscosity 
and  specific  heat  values  available  from  TPRC  or  other  reliable  sources. 
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